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SUMMARY 

This paper shows that it is within the capability of the 
DMAF language to build a set of vectors that can grow incremen- 
tally to be applied automatically and economically within a DMAP 
loop that serve to append sub-matrices that are generated within 
a loop to a core matrix. The method of constructing such vectors 
is explained. In an appendix are four DMAP packets for applying 
these techniques in different ways. 

PURPOSE 


The objective is to build a set of partitioning vectors 
that will grow incrementally to be applied automatically in a 
DMAP loop without analyst intervention for the purpose of append- 
ing the sub-matrices, formed in the loop, to the product matrix 
accumulated from preceding passes through the loop; i.e. 

T 

I ° 1 J 

* L 0 0 1 J T 

1^ 0 0 0 1 | T 


L 0 0 0 0 1 j T 
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DILEMMA 


The DMAP language in NASTRAN does not have a module that 
will generate vectors of arbitrary size. No matrix mathematical 
operations will increase the order of a matrix as demonstrated 
here . 

When two matrices are added or subtracted, the order of 
the result stays the same: 


"a b“ 

<t) 

'm n' 


'(a * m)(b ^ n>1 

- c d. 


-p q. 

2x2 

Jc + p) (d * q) j 


If two matrices are multiplied, the order of the result 
will not exceed the bounding order of the multiplicands: 


fa bl j"m n" 


"!am -i- bp) (an +■ bq)' 

UdJJ k 


_! cm + dp)(cn + dq)_ 


If a matrix is squared, the order of the result stays the 

same : 


“a £f 

n 

“(a^ + be) (ab + bd)' 

_c d_ 

2x2 

n 

_(ca + dc)(cb + d“i. 


Merging does produce a matrix of expanded order if one has 
at hand a partitioning vector, of the order desired, to monitor 
the merging. - This is tantamount to saying that you can expand 
the order of a matrix if you already have one around which is of 
the size that you want. A scheme is therefore needed which em- 
ploys a sequence of DMAP utility modules to provide an incre- 
mentally expanding partitioning vector. However, no single mod- 
ule available in NASTRAN will continually increase the order of 
a partitioning vector as desired here. If you let this notion 
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soak in for a time, then the inspiration percolates to the front 
of your hrain to say " start with soiuethincr bier and, extract the 
object of required order from it . " Our ideas at this point are 
still a long way away from the final algorithm, but this is the 
seminal idea. After a number of crude tries, the algorithm shewn 
here finally took shape. 

NOTATION 


A vector is a column matrix 



The transpose of a vector is a row matrix: 

,T 


tvr 


= j a l a 2 a 3 a 4 ' 


J 


To save space in the text, vectors will be written as transposes 
of row matrices: 



INITIAL CONSIDERATIONS 

The design of the DMAP ALTER packet will hinge on the 
range of certain parameters. 

1. The size of the core matrix to whioh sub-matrices are to be 
appended. (Number of starting columns or rows) 

2. The position at which sub-matrices are to be appended; i.e. 

Pre or Post. 

3. The size of the final matrix after all sub-matrices have been 
appended (Final number of columns or rows). 
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4. The number of columns or rows in the sub-matrix generated in 
each pass through the DMAP loop (Increment). 

SCHEME FOR APPENDING A SUB-MATRIX IN EACH PASS THROUGH A LOOP 

The nub of the idea is to operate with a series of vectors 
of order larger than that of the matrix to be built in the DMAP 
loop. The final MONITOR vector, that will be used for the actual 
merge operation, will be extracted from the SOURCE. The tech- 
nique which is used in the intermediate operations for managing 
the incremental growth of the MONITOR vector is where the real 
interest lies. This will be elaborated in detail, but as an 
introduction suffice it to say that there are two SOURCE vectors 
which are companions such that one has the increment at the tip 
of the vector and the other has the increment at the base of the 
vector. A BUBBLE vector is created whose length remains invari- 
ant, but whose increment moves inward from its starting position 
in a sea of zeroes. The BUBBLE vector leaves a trail of its pre- 
vious positions as an accumulation of the increments into a vec- 
tor called CLUSTER. The overall length of the CLUSTER vector 
remains invariant as the collection of increments increases. 
CLUSTER is used to partition the MONITOR from one of the SOURCE 
vectors . 


[SOURCE! .SHHBo CM0NIT0R3 

The steps in managing these vectors follow. For example, 
in the case of post-appending a single column in each pass, the 
SOURCE vector would have leading zeroes and a one in the last 
position, so a logical name would be BAS. 


|^0000 


0000000 


00000000000 


0001 J T 


Use this SOURCE vector as as a replenishable reservoir 
from which to take slices of increasing length to form the MONI- 
TOR vector as the product matrix grows from successive appending 
loops. 
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Several examples of these slices follow: 



Slice taken to create a MONITOR vector to merge a single vector 
to the right hand side of a core consisting of a single vector, 
resulting in a product matrix of 2 columns. 

00000001 j T 

Slice taken to create a MONITOR vector to merge a single vector 
to the right hand side of a 7-col core matrix, resulting in a 
product matrix of 8 columns* 

^ 000000000000001 J T 

Slice taken to create a MONITOR vector to merge a single vector 
to the right hand side of a 14-col core matrix, resulting in a 
product matrix of 15 columns. 

The CLUSTER vectors needed to do the slicing from the 
SOURCE would be of the same order as the SOURCE, but would have a 
cluster of increments equal to the number of columns in the pro- 
duct matrix after the end of a pass; i.e. the three CLUSTER 'vec- 
tors that would be used to slice the three MONITOR vectors shown 
above from the SOURCE would be respectively: 


000000 . 

0000000 . 

00000000000 . . 

0011 

J T 

000000 . 

ocooooo. . 

00000000000 

. .11111111 

'V 

r 

000000 . 


, 00000000000 111111111111111 

j T 

It is 

evident from the 

above three examples 

that the 

num- 


ber of ones in the CLUSTER vector grows by an increment after 
each pass (in this example the increment is one). The needs of 
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the algorithm will be met by saving only the latest CLUSTER vec- 
tor and devising some way of adding another increment in the next 
interior position. For example, in building the CLUSTER for the 
pass after the product has grown to order eight, one could add a 
vector to the "8" CLUSTER vector, which is of the same order as 
the "8" CLUSTER vector, but has zeroes everywhere except for a 
one in the 9th-f rom-the-end interior position. 

. .0000000 00000000 00011111111 j T 

* 

I 

+ I 

I 

0000000 ooooooooooo 00100000000 J T 

* 

Jr 

= Jr 

Jr 

0000000 00000000000 111111111 | T 

Well, here is a situation that can be used to evolve pass 
by pass. It is a matter of moving the "1" in the intermediate 
matrix to the next-interior-position . All well and good. But 
how does one do that? In a vulgar analogy, iust bite off a zero 
from the front and spit it onto the end. The tools are already 
at hand to do iust that. The two companion SOURCE vectors are of 
the right order. There is a ”1" on the front of one of the two 
companions to do the biting off of a zero from the front. The 
other of the pair has a "1" on the end to do the spitting of a 
zero onto the end. This intermediate vector is saved from pass 
to pass. Conceptually the "1" keeps travelling to the interior 
on successive passes. It is analogous to a bubble moving in a 
liquid-filled tube; thus the name BUBBLE vector. 
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Using this scheme, our needs will be met; 


1. The MONITOR vector grows by an increment on each pass to do 
the merging of an increment to the product. The MONITOR is dis- 
carded after each pass; and is regenerated on a succeeding pass. 

2. The CLUSTER vector augments its cluster of "l's" by an in- 
crement on each pass to do its job of slicing off a new MONITOR 
vector from the SOURCE vector. It is saved in its augmented form 
after each pass. 

3. The SUEBLE vector shifts its set of ones to the interior by 
an increment on each pass to do its job of adding an increment to 
the CLUSTER. It is saved after being shifted in each pass. 

4. The pair of SOURCE vectors remains unchanged during the whole 
operation. 

REQUIREMENTS 

The analyst is to supply simple initial bulk data and 
should be free of having to intervene in making entries in the 
DMAP code, once the author has incorporated the BUBBLE scheme 
into a DMAP loop. It should be easy for the author to use. It 
should be economical on storage space. 

ALGORITHM FOR APPENDING 

Now to build an algorithm to meet the stipulated require- 
ments, we will consider these requirements one by one. 

Specifications Simple initial bulk data. 

The user supplies 2 SOURCE vectors. Because they are 
sparse, it will take a minimum of 2 cards for each vector for a 
total minimum of four cards. This is simple. He knows how big 
"big" is from the nature of the problem that he is analyzing. He 
knows the size of the sub-matrix to be appended so he can set the 
vectors' increments. 
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Specification: Easy for the author. 

To illustrate the ease of concept, examples of the Ith 
step of a loop (1=4) will be used to demonstrate how simple it 
is to shift from pre or post appending or from column or tow 
appending. The operations can be characterized as CROP, SHIFT, 
SUM, HARVEST, AND APPLY. The first example shows pre-appending 
of columns and the second example shows post-appending of 
columns . 

Pre -appending 
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Specification: Economical on storage space 


The two source vectors will remain in storage for the 
whole time that the loop is in operation. Their length may be in 
the thousands, but the vectors are stored in packed form so the 
storage will consist of the header, the increment of l's, and the 
packing index for the position of nonzero entries. 

The CLUSTER will be the same length as the SOURCE vectors 
and will have a cluster of "l's" instead of a small increment of 
l's, but the content will be mostly zeroes in packed form; so its 
storage requirements will be low until its cluster grows to the 
final size of the PRODUCT matrix at the end of the loop. 

The BUBBLE vectors will be the same length as the SOURCE 
vectors and will have a single non-zero increment like the SOURCE 
vectors so their storage will be exactly the same as the SOURCE 
vectors . 


The MONITOR vectors will exist for the duration of a merge 
within a loop. Their brief storage requirements are exactly the 
same as the SOURCE vectors, because they will differ from the 
source only in the order number in the header. 

ALGORITHM FOR POST APPENDING A MATRIX 

Organize a pair of SOURCE vectors of large order; one with lead- 
ing zeroes before the increment of "l's"; and the other with 
trailing zeroes after the increment of "l's". Call these re- 
spectively BAS and TIP. 

BAS 1^00000 0 0 0 0 0 0 0 0 0 j I j T 

i The # of l's is 

governed by the sub-matrix order/pass; i.e. increment. 

TIP 1^00000 ....000000000 J T 

i The # of l's is governed by the increment/pass . 
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Find out the order of the core matrix prior to entering 
the loop. Use PARAML with operator = TRAILER and RECONO = 
trailer position # 1 corresponding to the number of columns in 

the core matrix. Call the parameter value that is returned by 
the name CORCOL for the number of core columns. But. if the sub- 
matrix is to be augmented by rows, set RECONO = trailer position 
#2 corresponding to the number of rows; then call the returned 
parameter by the name COROW for the number of core rows. The 
remainder of the explanation for this algorithm will be based 
upon the logic of post appending columns. 

PARAML C0REMTX//*TRAILER*/1 /CORCOL $ 

EQUIV COREMTX , PROD I /TRUE $ 

Move the bubble to that CORCOL position in unity incre- 
ments with these preparatory steps. 

Set up a small loop to prepare the vectors BUBLI and CLUSI 
ahead of the main loop. Companion source vectors (SOURCEl) with 
increments of unity called BAS1 and TIPI are supplied by the user 
with DMI cards in the bulk data. Make copies of BAS1 to act as 
starting configurations for CLUSI and BUBLI for molding into 
their loop configurations. 

COPY BAS1/CLUSI/0 $ For the case of post-appending. 

COPY BAS1/ BUBLI /0 $ 

Set the number of movements of the bubble that are needed 
to prepare BUBLI. 

PARAM / /*SUB* /SHIFT/ CORCOL /I $ 

LABEL BUBTOP $ 

FILE BUBLI= SAVE/ CLUSI = SAVE/ PRODI = SAVE $ 

Partition off a zero from the starting vector of BUBLI 
with TIPI. 
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C BUBLI} ==?!£!> { ■} 

(DUMMY) 


PARTN 


BUBLI. , TIPI /DUMMY,,,/ +7 $ 


Merge a zero onto the truncated stub of the BUBBLE vector 
with BAS1 


'DUMMY 


BAS1 


> CBUBLJJ 


MERGE 


DUMMY, , 


r 9 f 


BASl/BUBLJ / -t-7 $ 


The increment is now moved interior -ward in the BUBBLE vector. 
Combine the shifted BUBBLE vector with the starting vector of 
CLUSI . 


CCLUSI} + CBUBLJ} = CCLUSJ} 

ADD CLUSI .BUBLJ/CLUSJ/ $ 

Reset the internal loop names to the starting assignments in 
order to prepare for the next pass. 

SWITCH BUBLJ, BUBLI/ /-I $ 

SWITCH CLUSJ, CLUSI/ /-I $ 

Return to the starting label if the loop count (internally moni- 
tored by NASTRAN) is less than the prescribed value of SHIFT. 

REPT BUBTOP, SHIFT $ 

Now the positions of the ones in both BUBLI and CLUSI matches the 
size of the CORE matrix. 

Set up the loop for post appending an incremental sub-matrix onto 
the CORE matrix in each pass. 

LABEL CORTOP $ 
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Keep shifting the increment (which can now be different from 
unity depending on the need of the sub-matrix in the loop) in the 
BUBBLE vector once per' pass (with BAS and TIP instead of BAS1 and 
TIPI) now that it has been set by the preparatory packet. Start 
by truncating an increment of leading zeroes. 


CBUBLI1 ===i=> j 1 

Idummy ) 

PARTN BUBLI, , TIP/DUMMY , , , /+7 $ 

Merge on an increment of trailing zeroes to push the bubble into 
the interior. 


/DUMMY 
1 0 
MERGE 

Combine the 
"CLUSI". 


=i=§=> CBUBLJ3 

DUMMY, , , , 
shifted BUBBLE 


,BAS/BUBLJ/+7 $ 
vector with the 


CLUSTER vector 


CCLUSI} + CBUBLJ1 = CCLUSJJ 


ADD CLUS I , BUBL J / CLUS J / $ 

Partition off the MONITOR vector to the size of the PRODUCT 
matrix being built in this pass. 

cbasj .SWH.> I-!-} 

ImntriJ 

PARTN BAS, , CLUS J / , MNTRI , , /+7 $ 

The desired merging vector has been produced and is ready to be 
applied to the job of joining the current sub-matrix to the cur- 
rent PRODUCT matrix. 
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$$MOCK FUNCT'L MODULE "MTXOPN" PRODUCES THE SUB-MATRIX "SUBMTX" 
MTXOPN - DB IN /SUBMTX/ $ 

$$ PRODUCT MATRIX "PRODI" IS AVAILABLE FROM A PREVIOUS PASS. 

CPRODI : SUBMTXD ==^E£i> CPRODJD 

MERGE PRODI, , SUBMTX, ,MNTRI , /PRODJ/ +7 $ 

Prepare for the next pass by resetting the matrix names to their 
initial assignments. 

SWITCH PRODJ, PRODI/ /-I $ 

SWITCH CLUSJ , CLUSI / /-I $ 

SWITCH BUBLJ,BUBLI/ /-I $ 

Apply a test to see how far the merging has proceeded, with re- 
spect to the requirements of the matrix operations. Set the 
value of a threshold to negative when the loop is done so that a 
conditional jump takes the OSCAR outside the loop. 


$$MOCK FUNCT'L MODULE "TESTOPN" PROVIDES A TEST LOGIC FOR LOOPS. 
TESTOPN VARY, , / DBOUT / THRESHLD S 

COND OUTS IDE, THRESHLD $ 


Loop back to the top of the main loop so long as the value of 
THRESHLD remains positive. When it turns negative the condi- 
tional jump takes the sequence of operations beyond the loop. 


$$ 

$$ 


REPT 

LABEL 


CORTOP,## $ ## VALUE IS LARGER THAN THAT 
EXPECTED FOR THRESHOLD 

OUTSIDE $ PRODI IS PASSED OUTSIDE THE LOOP WITH 
SUB -MATRICES FROM ALL LOOPS MERGED ONTO IT 
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This completes the explanation for the algorithm following 

the example of”' post-appending any sized sub-matrix by columns. 

In the appendix, this algorithm is coded for the four permuta- 
tions of : 

1. Post-append by columns. 

2. Pre -append by columns. 

3. Post append by rows. 

4. Pre -append by rows. 

CONCLUSION 

It has been shown that it is within the capability of the 
DMAP language to build a set of vectors that can grow incremen- 
tally to be applied automatically and economically within a DMAP 
loop to serve to append sub-matrices that are generated within a 
loop to a product matrix. 


131 



APPENDIX FOR FOUR PATHS OF THE BUBBLE ALGORITHM 


$POST APPEND COLUMNS. 


PARAML 

COREMTX/ /^TRAILER*/ 1 /CORCOL 

At 

COPY 

BAS1 /CLUSI/ 0 $ 


COPY 

EAS1 /BUBLI/ 0 $ 


PARAM 

/ / *SUB* / SHIFT / CORCOL / 1 $ 


LABEL 

BUBTOP $ 


FILE 

BUBLI = SAVE/ CLUSI = SAVE/ PRODI 

PARTN 

BUBLI, , TIPI /DUMMY, , , /+7 

$ 

MERGE 

DUMMY, , , , , BAS1 /BUBLJ/+7 

$ 

ADD 

CLUSI , BUBLJ /CLUS J / $ 


SWITCH 

BUBLJ, BUBLI/ /-I $ 


SWITCH 

CLUSJ, CLUSI/ /-I $ 


REPT 

BUBTOP, SHIFT $ 


LABEL 

CORTOP $ 


PARTN 

BUBLI, ,TIP /DUMMY, , , .'+7 

$ 

MERGE 

DUMMY, , , , , BAS /BUBLJ/ +7 

$ 

ADD 

CLUS I, BUBLJ /CLUSJ/ $ 


PARTN 

BAS , , CLUSJ/ , MNTR I , , /+7 

$ 

MTXOPN 

DBIN/MODE/ $ 


MERGE 

PRODI, ,MODE , ,MNTRI , /PRODJ/ 

SWITCH 

PRODJ, PRODI/ /-I $ 


SWITCH 

CLUSJ, CLUSI/ /-I $ 


SWITCH 

BUBLJ, BUBLI/ /-I $ 


TSTOPN 

VARY, , / DBOUT / THRESHLD $ 


COND 

OUTS IDE, THRESHLD $ 


REPT 

CORTOP, NMBR $ 


LABEL 

OUTSIDE $ 



SAVE $ 
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$PRE APPEND COLUMNS 


PARAML 

COPY 

COPY 

PARAM 

LABEL 

FILE 

PARTN 

MERGE 

ADD 

SWITCH 

SWITCH 

REPT 

LABEL 

PARTN 

MERGE 

ADD 

PARTN 

MTXOPN 

MERGE 

SWITCH 

SWITCH 

SWITCH 

TSTOPN 

COND 

REPT 

LABEL 


COREMTX/ / ^TRAILER* / 1 ,'CORCOL $ 

TIPI /CLUSI/ 0 $ 

TIPI /BUBLI/ 0 $ 

/ /*SUB*/SHIFT/CORCOL/ 1 $ 

BUBTOP $ 

BUBLI = SAVE/ CLUSI = SAVE/ PRODI = SAVE $ 
BUBLI, , BAS1 /DUMMY, , , /+7 $ 

DUMMY, , , , .TIPI /BUBLJ/ +7 $ 

CLUSI ,BUBLJ /CLUSJ / $ 

BUBLJ, BUBLI/ /-I $ 

CLUSJ, CLUSI/ /-I $ 

BUBTOP, SHIFT $ 

CORTOP $ 

BUBLI, ,BAS /DUMMY, , , /+7 $ 

DUMMY, , , , ,TIP /BUBLJ/ +7 $ 

CLUSI, BUBLJ /CLUSJ/ $ 

TIP , .CLUSJ/ ,MNTRI , , /+7 $ 
DBIN/MODE/ $ 

PRODI, .MODE, ,MNTRI, /PR0DJ/+7 $ 

PRODJ, PRODI/ /-I $ 

CLUSJ, CLUSI/ /-I $ 

BUBLJ, BUBLI/ /-I $ 

VARY, , / DBOUT / THRESHLD $ 

OUTS IDE, THRESHLD $ 

CORTOP, NMBR $ 

OUTSIDE $ 
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$POST APPEND ROWS 


PARAML 

COPY 

COPY 

PARAM 

LABEL 

FILE 

PARTN 

MERGE 

ADD 

SWITCH 

SWITCH 

REPT 

LABEL 

PARTN 

MERGE 

ADD 

PARTN 

MTXOPN 

MERGE 

SWITCH 

SWITCH 

SWITCH 

TSTOFN 

COND 

REPT 

LABEL 


COREMTX/ /^TRAILER*/ 2/COROW $ 

BAS1 /CLUSI/ 0 $ 

BAS1 /BUBLI/ 0 $ 

/ / *SUB*/ SHIFT/ COROW/ 1 $ 

BUBTOP $ 

BUBLI = SAVE/ CLUSI = SAVE/ PRODI = SAVE $ 
BUBLI, , TIPI /DUMMY, , , ,'+7 $ 

DUMMY, , , , , BAS1 /BUBLJ/+7 $ 

CLUS I , BUBL J / CLUS J / $ 

BUBLJ, BUBLI/ /-I $ 

CLUS J, CLUSI/ /-I $ 

BUBTOP, SHIFT $ 

CORTOP $ 

BUBLI, ,TIP /DUMMY, , „ / f 7 $ 

DUMMY, , , , , BAS /BUBLJ/ +7 $ 

CLUS I, BUBLJ/ CLUS J/ $ 

BAS , , CLUSJ / , MNTR I , , ,'+7 $ 

DB IN /MODE/ $ 

PRODI, MODE , , , , MNTR I /FRODJ/+7 $ 

PROD J, PRODI/ /-l $ 

CLUSJ, CLUSI/ /-I $ 

EUELJ ,EUELI/ /-I $ 

VARY, , / DBOUT / THRE3HLD $ 

OUTS IDE, THRESHLD $ 

CORTOP, NMBR $ 

OUTSIDE $ 
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$PRE APPEND ROWS 


PARAML 

COREMTX/ / ^TRAILER* / 2 / COROW $ 


COPY 

TIPI /CLUSI/ 0 $ 


COPY 

TIPI /BUBLI/ 0 $ 


PARAM 

//*SUB*/ SHIFT/COROW/ 1 $ 


LABEL 

BUBTOP $ 


FILE 

BUBLI = SAVE/ CLUSI = SAVE/ PROD I 

PARTN 

BUBLI, , BAS1 /DUMMY, , , /+7 

$ 

MERGE 

DUMMY, , , , ,TIP1 /BUBLJ/+7 

$ 

ADD 

CLUSI , BUBL J / CLUS J / $ 


SWITCH 

BUBLJ, BUBLI/ /-I $ 


SWITCH 

CLUSJ, CLUSI/ /-I $ 


REPT 

BUBTOP, SHIFT $ 


LABEL 

CORTOP $ 


PARTN 

BUBLI, , BAS /DUMMY, , , /+7 

$ 

MERGE 

DUMMY, , , , ,TIP /BUBLJ/ +7 

$ 

ADD 

CLUS I, BUBLJ /CLUSJ/ $ 


PARTN 

TIP , , CLUSJ/ ,MNTRI , , /+7 

$ 

MTXOPN 

DBIN/MODE/ $ 


MERGE 

PRODI ,MODE, , , ,MNTRI /PRODJ/ 

SWITCH 

PRODJ, PRODI/ /-I $ 


SWITCH 

CLUSJ, CLUSI/ /-I $ 


SWITCH 

BUBLJ, BUBLI/ /-I $ 


TSTOPN 

VARY, , / DBOUT / THRESHLD $ 


COND 

OUTS IDE, THRESHLD $ 


REPT 

CORTOP, NMBR $ 


LABEL 

OUTSIDE $ 



SAVE $ 
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